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We introduce a measure Q of bipartite quantum correlations for arbitrary two-qubit states, expressed as a state- 
independent function of the density matrix elements. The amount of quantum correlations can be quantified 
experimentally by measuring the expectation value of a small set of observables on up to four copies of the state, 
without the need for a full tomography. We extend the measure to 2 x d systems, providing its explicit form 
in terms of observables and applying it to the relevant class of multiqubit states employed in the deterministic 
quantum computation with one quantum bit model. The number of required measurements to determine Q in 
our scheme does not increase with d. Our results provide an experimentally friendly framework to estimate 
quantitatively the degree of general quantum correlations in composite systems. 
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Quantum entanglement is one of the most fundamental 
consequences of the superposition principle and undoubtedly 
plays a key role in designing faster-than-classical algorithms, 
teleportation protocols and super-dense coding (TJ. However, 
it may not be the ultimate resource behind the power of quan- 
tum computation [2|. It has been recently found that, even 
with no entanglement, some mixed-state based schemes [such 
as the so-called deterministic quan- turn computation with one 
quantum bit (DQC1) [3]] allow an improvement of perfor- 
mance in computing tasks [4], and, more generally, separa- 
ble states possess genuinely quantum correlations [5| (QCs), 
captured e.g. by the quantum discord J6] [7J , which cannot be 
described within a classical scenario II6 UTTI . In general, QCs 
in a state pab can be defined as the minimum amount of to- 
tal correlations (measured, e.g., by the mutual information) 
between Alice and Bob, that are destroyed by a local mea- 
surement on one or both subsystems 1161481 [lOl [T2l4T5ll . For 
pure states, QCs coincide with entanglement [;6i|. For mixed 
states, even if some operational interpretations have been pro- 
posed lfTTl[T5|jT8l , basic technical issues still prevent us from 
reaching a full comprehension of their nature. Indeed, theo- 
retical evaluation and experimental detection of QCs both rep- 
resent hard challenges: any attempt to determine the QCs in a 
given state pab is hindered by the difficulty of solving an op- 
timisation to determine the least disturbing measurement for 
that state, thus requiring the full knowledge of it. Recently, 
some non-tomographic detection schemes for witnessing non- 
vanishing QCs by measuring just one observable lfl4l [T9l 1201 
have been proposed and implemented. However, noting that 
all states possess nonzero QCs but a null-measure set ETI . 
the most worthwhile question becomes that of evaluating, by 
a proper measure, the actual amount of QCs encoded in a state. 
Only then, quantitative connections can be drawn between the 
QCs content and the performance of some quantum protocol 
using them as resource HI l22l l23l . 

In this Letter, we show that QCs in a general two-qubit state 
p can be reliably quantified without any explicit optimisation 
and with no need to know the full shape of the state. We define 
a QCs measure Q which is a state-independent function of the 
density matrix elements. In particular, Q can be expressed in 



terms of the expectation values of a set of nine observables 
{Oj}. Consequently, such a function could be evaluated by 
designing simple quantum circuits simulating the measure- 
ments of \{Oi)} [24-31]. However, following the alternative 
approach of Refs. [32-34], we further show that the quan- 
tity Q can be even less demandingly measured by perform- 
ing seven local projections on up to four copies of the state 
p . Then, we extend our measure to capture bipartite QCs in 
states of 2 x d dimensional systems, finding that seven pro- 
jective measurements are always sufficient to experimentally 
determine Q, i.e. the number of measurements required is in- 
dependent of d. Specifically, we use this construction to ob- 
tain a quantitative estimate of QCs in a recent experimental 
implementation [35| of the DQC1 model [3| with four qubits. 

A number of conceptually different measures of general 
QCs have been recently proposed HHE] [TUl IT214T51 . In the 
following, we consider a two-qubit state p = pab and adopt 
a geometric perspective, quantifying the QCs in terms of the 
minimum distance of p from the set Q. of classical-quantum 
states. The states ^ € filling such set are left unperturbed 
by at least one choice of projective measurement on Alice, and 
take the form IfTTll x — Pi\i)(i\ ® PiB, where p-, are probabil- 
ities, {\i}} is an orthonormal vector set and p,B is the marginal 
density matrix of Bob. Adopting the Hilbert-Schmidt norm 
||M|| 2 = VTr^MMt), one obtains a QCs measure known as 
"geometric discord", introduced in [14], operationally inter- 
preted in lfT71l23l . and defined as 

D G (p) = 2 rnin \\p (D 

where we add a normalisation factor 2. The geometric dis- 
cord enjoys a closed expression for two-qubit states. First, 
one needs to express the state in the Bloch basis, p = 

\ Ey=o ® o-j = i(i4 + 2 L w ® h + E;=i yjh ® o-j + 

2y=i tijo-i ® o-j\ where R u = Tr[p(o-,- ® cr ; )], cr = I 2 , cr ; 
(z = 1,2,3) are the Pauli matrices, x = ]x,],y = {y,} repre- 
sent the three-dimensional Bloch column vectors associated to 
A, B, and fy are the elements of the correlation matrix t. Then, 
following 01, we have: D G {p) = ±(||^| 2 + \\t\\\ - 4£ max ) = 
2Tr[5] - 2fc max , with fc max being the largest eigenvalue of the 
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FIG. 1: Geometric discord D G versus its tight lower bound Q for 
3 • 10 4 random two-qubit states. The plotted quantities are dimen- 
sionless. 

matrix S = ^(xx J + ff T ). We now provide an explicit expres- 
sion for £ raax . The characteristic equation of the matrix S is 
cubic and can be solved analytically l36ll . Being constrained 
to real solutions only, we write the eigenvalues {£,} of S as 

_ Tr[S] V6Tr[S 2 ] -2Tr[S] 2 (O + aA 



where {a,} = {0,2^, 4n} and 9 = arccos j(2Tr[S] 3 - 
9Tr[S]Tr[S 2 ] + 9Tr[S 3 ]) V2/(3Tr[S 2 ] - Tr[S] 2 ) 3 }. Since 9 is 
an arccosine, we have < 9/3 < n/3 and the maximum of 
cos ( ^- L is reached for or; = a\ — 0. Hence, k max = max{£,} 
k , and the geometric discord for an arbitrary two-qubit state p 
assumes the form of a state-independent function of its entries 
(fij), that is 

D G (p)=2(Tr[S]-h). (3) 

However, we aim to define a simpler, and more accessible ex- 
perimentally, QCs quantifier. By replacing 9 with in Eq. Q, 
we obtain a meaningful and remarkably tight lower bound 
Q < D c (see Fig.[T| to the geometric discord, given by 

Q(p) = | (2Tr[S ] - V6Tr[S 2 ]-2Tr[S] 2 ) . (4) 

The quantity Q is still a state-independent expression of the 
entries of p, but a rather easier one to manage than Dq, and 
can be regarded as a bona fide measure of QCs in its own 
right. Indeed, it is non-negative by definition, it is faithful 
(i.e., vanishes only on classical-quantum states x) an d coin- 
cides with D c for pure states. The latter two properties can 
be proven as follows. Faithfulness is equivalent to showing 
that Q = <==> Dq = 0; the condition for vanishing Q 
is Tr[5] 2 = Tr[5 2 ], and by the Cayley-Hamilton theorem this 
implies Tr[S] 3 = Tr[5 3 ], i.e., Dq = 0. The equality between 
Dq and Q for pure states follows from the fact that, writing 
a bipartite pure state (where subsystem A is a qubit) in the 
Schmidt decomposition, we have 9-0. We also find that Q 
provides a non-trivial upper bound on an entanglement mea- 
sure, specifically the squared negativity TV 2 [1|. In fact, the 
chain Dq > Q > N 2 holds for arbitrary two-qubit states, with 
all inequalities saturated on pure states [39]. Finally, let us 
mention that a simple upper bound on Dq can be obtained as 
well from Eq. 0, D c < 4Tr[5]/3. 



From now on, we adopt Q as a rightful QCs quantifier 
for two qubits, endowed with the advantage of requiring nei- 
ther theoretical optimisations nor experimental state tomog- 
raphy for its evaluation. Specifically, the task of providing 
a recipe for measuring Q reduces to write Ti[S] and Tr[5 2 ] 
as functions of suitable observables, and is accomplishable 
as follows. Defining the matrices X = xx 1 , T = tt T , we 
have Tr[S] = (Ti[X] + Tr[T])/4 = Tr[p 2 ] - Tr[p|]/2, and 
Tr[5 2 ] = yg(Tr[X 2 ] + Tr[r 2 ] + 2Tr£XT]). After some algebra 
we obtain 

Tr[S 2 ]= | (-2 - 8Tr[p 4 ] + 8Tr[p 3 ] + 6Tr[p 2 ] 2 (5) 

- 2Tr[p 2 ](5 + Tr[p 2 ]) - 2Tr[p 2 ] 2 + 10Tr[p 2 ] 

- Tr[p 2 ] 2 + 12Tr[p 2 ] - 6Tr[p 2 ]Tr[p|] 

+ 4Tr[p(I 2 ®p B )p(I 2 ®p fl )] - 24Tr[p(p A ®p B )] 
+ 8Tr[p(p A 9 l 2 )p(p A ® I 2 )] + 8Tr[p 2 (p A ® p B )]). 

By substituting in Eq. Q, Q takes the form of a function 
of fourth-order polynomials of (p,y), in particular it is writ- 
ten in terms of traces of matrices powers. Now, given a 
general density matrix p, it holds |[24H27l [30l |3~TH Tr[p*] = 
Tr[VV®*], where V k is the shift operator, VVi^...ft) = 
l^ifc^i . . - <pk-\)- Also, for two unknown states pi,p 2 , it has 
been proven [27 1 that Tr[V 2 pi ®p 2 ] = Tr[pip 2 ]. More gener- 
ally, we have 1371 Tr[pip 2 . . .p,t] = Tr[V^p! ® p 2 ® . . . ® p^]. 

We can exploit these results and follow the approach pro- 
posed in ll27l for estimating the expectation values of the ap- 
propriate unitary operators {0;}? =] to associate with each of 
the nine independent factors in Eq. |5]l (which include those 
appearing in Ti[S] as well). They can all be expressed as shift 
operators 0-, = V k on a number k (k < 4) of copies of the 
global and/or marginal density matrices and their overlaps, 
depending on each particular term in Eq. |5]). The circuit to 
be implemented, which includes an ancillary meter qubit, is 
depicted in Figj2ja). For each 0„ we build an interferometer 
modified by inserting a controlled- (9, gate: defining the visi- 
bility v, we obtain in general Tr[(9,pi®p 2 ®. . .®pk\ = v. Hence, 
QCs can be measured quantitatively from the expectation val- 
ues of the nine operators {(9,) only — as opposed to 15 ob- 
servables required for complete state tomography — obtained 
in laboratory via readouts on the ancillary qubit. 

We wish now to provide an alternative scheme for the exact 
measurement of Q that further reduces the required resources 
for its implementation, and appears even more experimentally 
friendly. This is done by rephrasing the detection scheme in 
terms of local (with respect to the Alice-Bob split) projectors 
on multiple (up to 4) copies p m of the same state p ll28l l29l 
I32H341 . We observe that 

Tr[S] =Tr[p 2 ]-Tr[p|]/2. (6) 

It is known El thatTr[p 2 ] = Tr[V 2 p® 2 ] = Tr[(P + -f )p® 2 ] = 
1 - 2Tr[P~p® 2 ], where V 2 is the swap operator and are 
the projectors on the symmetric/antisymmetric subspaces. 
Named A,- (Bj) the subsystems controlled by Alice (Bob) in 
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the z'-th O'-th) copy of the bipartite state p, we have then 
1 



Tr[S] 



ZTrFcUxW'H + T ^b 2 Pb 1 . (7) 



where for two qubits P B ^ = I^^X^J, |^ l2?2 > = 
-L(|01>-|10», while ^ lBl)(A2B2) = |(3I 16 -2;^ 4) <8>^ 4) ), 
where with cr (d) we indicate the generalised (and normalised) 
Gell-Mann matrices for dimension d. Alternatively, we can 
exploit the very recent results of Ref. P4ll and write 



Tr[S] 



4ci - 2c 2 - c 3 + 



d = Tr[(P AiA2 
ci = Tr[(P AiA2 
c 3 = Tr[(I AlA2 ®P„ & )p 82 ] 



I BlB2 )p 82 ] ; 



(8) 



In 11291 a method to measure the purity of a quantum state, 
which is all that we need, is presented and demonstrated by 
means of the implementation of an all-optical set up. Ref. l38l 
presents a more comprehensive detection scheme for projec- 
tive measurements; see also 11401 for a very recent alternative 
method. To sum up, in this framework we need three mea- 
surements of two-qubit projectors [Eq. ([8]l] — or two measure- 
ments, one on two qubits, the other (non-local with respect to 
the Alice-Bob split) on four qubits [Eq. ([7]i] — and two copies 
of the state, to measure Tr[S ]. 

The detection of Tr[5 2 ] can also be recast in terms of local 
projections. Following Ref. [34j, we obtain 



Tr[S 2 ] 



C4 



16c 4 + 8(c 7 - c 5 - 2c 6 ) + c\ + \c\ - c 3 -2c 2 + \ 



TrtC^AtAi w 1 A 2 A 

Tr[(P AiA4 ®I A2A3 
Tr[(I AlAi ®P A , 



)P m ] ; 



9P Bl B 2 ®h iBi )p® 4 ]; 
>P Bi B 2 ®h,B 4 )p m ]. 



(9) 



Compared to the measurement of Tr[5], here we have again 
projectors on pairs of qubits, however they need to be imple- 
mented on four copies of the state p. As Eq. Q shows, we 
can evaluate the value of Tr[5 2 ] by measuring four such inde- 
pendent projectors in laboratory. Therefore, in the most eco- 
nomical scheme devised here, the full quantitative detection of 
bipartite QCs in an arbitrary two-qubit state p as measured 
by Q demands six or seven projective measurements on ( up 
to) four copies of the state p. Notice for comparison that to 
measure the geometric discord D c exactly, Eq. <|3j, one would 
need 1 1 projective measurements on up to six copies of the 
state |3~4*1 . On the other hand, at a qualitative level, a single 
observable witness suffices to reveal whether Q (or the dis- 
cord) is zero or not lfl4l [191 l20l l35l . 

We now extend our measure to higher-dimensional sys- 
tems, in particular to 2 X d systems, which include the prac- 
tically relevant case of one qubit (A) versus a register (B) of 
n qubits. The geometric discord for an arbitrary state p of 
a 2 x d system has been derived in RTl . and has the same 
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FIG. 2: (a) Quantum circuit estimating Tr[0,pi ®P2 ® . . . ® Pk\ = v. 
Two Hadamard gates H are applied to an ancilla, followed by a mea- 
surement in the computational basis. The interferometer is modified 
by inserting a controlled-O, gate acting on the overlap of states, 
(b) DQC 1 model with a register of three maximally mixed qubits and 
an ancillary qubit in a state of polarisation \i. Expectation values of 
<Ti , <r 2 on the ancilla return the real and imaginary part of Tr[(/]/8. 



expression as Eq. ([3]), just amended with the following gen- 
eralisations: Tr[5] = 5b(Tr[X] + Tr[!T]), x { = Tr[p(o-i ® I d )], 
and tjj = Tr[p(cr ; ® Ty)], where, for example, we can assume 

{r ,} = {<T (d) } as af-dimensional basis for Bob's subsystem. We 

J j 

can repeat the steps done for 2x2 systems and obtain a state- 
independent form for Dc and for the lower bound Q as well, 
since the subsystem A is still a qubit and therefore the charac- 
teristic equation of S remains a cubic. Thus, for 2 x d states, 
the task we face is again to express Tr[5],Tr[5' 2 ] in terms of 
observables. 

The most practical way to proceed is to consider the scheme 
in terms of local projectors. In this respect, it is straight- 
forward to verify that Eqs. Q ^ still hold for 2 x d sys- 
tems: their expression can be written in exactly the same form 
as for the two-qubit case, provided we generalise the swap 
and the projectors P to arbitrary dimension d as follows, 
V 2 = L d (\p + E; T f ® Tt) and P" St = ±{{d - Vfo - Tj ® T t ), 
where S j,Sk denote two <i-dimensional systems, and the t,'s 
reduce to Pauli matrices in dimension d — 2 (e.g., when we 
want to calculate Trfp^] in the 2x2 case). This observation, 
combined with the previous analysis, allows us to conclude 
that, even/or arbitrary states pof2xd dimensional systems, 
we just need six or seven projective measurements on up to 
four copies of the state p to quantify bipartite QCs between 
the qubit and the remaining qudit system. The number of mea- 
surement settings thus does not increase with d, which demon- 
strates the efficiency and scalability of our scheme. Clearly, 
the optical implementation of projectors of the type PbiBj, 
i.e. multi-qubit projectors, is more complicated than the two- 
qubit case, see e.g. 02). However, the method demonstrated 
in 1 38 1 can be extended to arbitrary dimensions without dra- 
matically increasing the complexity of the experimental set- 
ting (as claimed by the authors in the last section of Ref. Il38ll ). 
More precisely, the number of optical elements required to 
implement each projector (basically interferometers) should 
increase poly normally — namely, linearly — with d l38ll43l . in 
stark contrast with a complete quantum state tomography for 
which the required resources scale exponentially B4l . Note 
also that our scheme for 2 x d systems is completely general 
and no prior knowledge of the form of the state is required; 
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FIG. 3: Bipartite quantum correlations as measured by Q (red con- 
tinuous line) and D a (black dashed line) for the output state of the 
DQC1 model as implemented in Ref. [35] with four qubits, plotted 
as functions of the initial polarisation /j of the ancilla qubit. All the 
plotted quantities are dimensionless. 

it only relies on the implicit assumption that the subsystem A 
has dimension 2, i.e., it is indeed a qubit. This assumption can 
be verified in laboratory a priori e.g. by measuring suitable 
Hilbert space dimension witnesses 14511 . or, possibly, with to- 
mography on the marginal state of subsystem A, which only 
consumes a fixed, small amount of extra resources. 

We consider as an example the four-qubit implementation 
of the DQC1 model for quantum computation [3 1. This algo- 
rithm estimates the trace of a normalised unitary matrix U. We 
consider the instance recently implemented experimentally in 
11351 (where only a discord witness rather than a quantitative 
estimate was measured), where U = (a, a,b,\,a,b,\, 1), with 
a = -(e - ' 3 ^ 5 ) 4 , b = (e~' 3lr / 5 ) 8 . Such a specific gate is used for 
the approximation of Jones polynomials |46|. The first qubit 
(the ancilla) is initially in a state of polarisation ju, while the re- 
maining qubits are maximally mixed. Referring to the scheme 
of Fig. [2jb), the final state of the system before readout is 
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Is 



(10) 



We calculate the bipartite QCs Q(p ou i), measurable in labo- 
ratory according to the scheme detailed above, between the 
ancilla and the residual three-qubit system, and we compare 
them with the geometric discord Dq, while entanglement is 
always zero across this bipartition. The plots in Fig. [3]reveal 
that Q is in good agreement with Dq, being a mo no tonic func- 
tion of the polarisation p. of the ancilla, hence showcasing its 
reliability as a QCs quantifier (4). 

In conclusion, we presented a scheme to quantify theoret- 
ically and experimentally general bipartite QCs for arbitrary 
two-qubit and qubit-qudit states. We introduced a measure 
Q that is a state-independent function of polynomials of the 
density matrix elements, and can be measured by implement- 
ing a restricted number of quantum circuits, or alternatively a 
restricted number of local projections, on up to four copies 
of the state, which appears in reach of current technology 
Il33ll38l 43). We used our measure to evaluate quantitatively 
the degree of QCs created in a recent experimental implemen- 
tation [35 1 of the DQC1 model with four qubits |3). 

Providing experimentally friendly recipes for the measure 
of QCs in n-partite realisations of quantum information pro- 



tocols is key to clarifying their usefulness for the performance 
of such practical tasks |2). In this respect, much attention is 
being devoted to the QCs dynamics in open quantum systems 
R71 -I501 and, independently, to characterising the transition 
from Markovian to non-Markovian regimes |5T}[55). Non- 
Markovianity can be witnessed by monitoring entanglement 
between one subsystem, coupled to the environment, and an- 
other clean subsystem 11541 . One might imagine that more 
general QCs could be somehow more sensitive to the proper- 
ties of dynamical maps. However, it is known that even local 
Markovian channels (as well as Markovian common environ- 
ments [56]) can induce an increase of discord-like QCs in a 
composite system l57l . Therefore, the question needs to be 
formulated properly and with care, and demands a dedicated 
analysis which is beyond the scope of this work |58|. 

We hope our Letter may contribute to render the general 
quantumness of correlations a more accessible (theoretically 
and experimentally) concept in the study of complex quantum 
systems. 
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